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SYNOPSIS 


A Lheois entitled Qualitative Properties of Solutions of 
Volterra Integro~Differential Equations" is sutraitted in 
partial fulfilment of the requirement for the Ph.D. degree 
by Pyda Srinivas to the Department of Mathematics, Indian 
Institute of Technology, Kanpur* 

Quite often linear and nonlinear Volterra integro— 
differential equations arise in diverse areas such as 
physical, biological, social and engineering sciences* The 
way such equations arise and the properties of solutions 
of these equations to various physical phenomena have been 
investigated by many scientists. However, until quite 
recently, the attempts have not been made to investigate the 
qualitative properties such as boundedness, oscillatory and 
stability of these equations. The object of this thesis is 
to study the positivity, boundedness, eventual stability 
and L^-stability properties of solutions of Volterra integro- 
differential equations. Dichotomy properties of linear 
equations also find a place in this thesis. 

Friedman,Komlenko , Levin, Miller and Padmavalli among 
others have studied the positivity and boundedness of solutions 
of Volterra integral equations and many interesting results 
have been reported. In the case of integral equations of the 


form 
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t 

x(t) = f(t) - I a(t,s) x(s)ds^ O < t. < w (1) 

o 

one can assert by virtue of Netmann series that all the solutions 
of (1) are positive in their interval of existence, if a(t,s) 

IS negative and f(t) is positive. However, the same is not 
true in general for integro— differential equations of the 
form 

t 

x'(t) = f(t) - / a(t,s)x(s) ds, x(0) = x^^ > O (2) 

o 

J>articularly -vdien a(t, s) is positive. Sufficient conditions 
for positivity and boundedness of solutions of (2) are 
investigated in the case >iiere a(t/s) is positive, CXir 
approach here involves the construction of a scalar function 
which in certain cases turns out to be the usual resolvent 
kernel. Examples and remarks are provided at appropriate 
places to indicate the advantage of this approach and 
illustrate the generality of conditions. 

Eventual uniform asymptotic stability is a natural 
generalisation of uniform asymptotic stability of a given 
differential equation in v^ich it is not assumed that the 
zero function is a solution of it. This situation would in 
general happen for integro— differential equations of the foxan 

t 

x^(t) = f(t) + / a{ t,s)G( s,x(s) )ds (3) 

o 

When either f(t) f O or G(t/0) ^ O, In the later part of 
the thesis, we shall investigate sufficient conditions for 
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eventual uniform asymptotic stability of the set {0| for (3) 
when G{s,x) is linear in x* The results thus obtained are 
used to study the eventual properties of solutions of the 
nonlinear equation (3) under more general class of non-linear 
functions G. 

Liapunov second method involving the construction of a 
scalar function is one of the powerful tools in the qualitative 
study of ordinary differential equations and Volterra integro- 
differential equations. Concerning the application of this 
method for functional differaatial equations^ essentially two 
approaches have emerged in recent years such as the 
construction of a Liapunov (i) function and Cii) functional# 

The Liapunov^Razumikhin technique is basically of the former 
type and this has been adapted to discuss the L -stability 
properties of solutions of Volterra integro-differential 
equations. 

The study of linear analysis plays a vital role to 
understand the complexity of the corresponding nonlinear 
analysis. Most of the authors in their recent studies on linear 
Volterra integro-differential equations either by Liapunov 
second method or from perturbation analysis^ have assxamed 
that the coefficient matrix is stable. Notable exceptions 
that have dispensed with this condition on the coefficient 
matrix have been the investigations of Levin, GrosOTtan and 
Miller, Grimmer and Seifeirt, and Burton among others# Motivated 



iv 


from the intearestrng nature of the problem^ an attCT^st has 
been made in this thesis to study exponential stability 
properties of solutions of Volterra integro— differential 
systems when the coefficient matrix is not necessarily stable. 
Our approach here has been to develop an equivalent equation 
which involves an arbitrary function. The proper choice of 
this function would pave a way for the new coefficient matrix 
(corresponding to equivalent Volterra integro— differential 
equation) to be stable. The stability analysis of linear 
autonomous systems basically characterises the sign of the 
real parts of the characteristic roots of the coefficient 
matrix. However, this approaches does not carry over to 
non-autonomous systems. Thus the study of dichotomies of 
solutions of linear non-autonomous systems assximes a 
prominent role. In the last chapter of this thesis the 
dichotomy properties of linear Volterra integro— differential 
systems are discussed. 



CHAPTER - 1 


INTRODUCTION 

1 .1 Historical Notes 

Ever since the beginning of calculus, many physical 
phenomena have been modelled as equations involving derivatives 
or integrals. Naturally as a follow-up of any modelling, it 
is necessary that one should know about the solutions of 
these equations. During the initial decades of its existsice, 
this theory consisted only Of isolated methods of solving 
certain type of linear equations. The problem of existence 
of a solution and its explicit representation was posed 
already in the eighteenth century itself. Numerous 
investigations revealed that the e3q)licit representation of 
the solutions is an extremely rare phenomenon and in many 
cases (particularly for linear equations with variable 
coefficients and nonlinear equations), it is rather difficult 
to obtain the same in the e3q>licit form. Further, the methods 
of iiumerical integration of these equations did not open the 
horizons for the development of the general theory since 
these methods yield only one particxilar solution that too 
pn a finite interval. Applications — in particular to the 
problems of celestial mechanics, bio-mechanics, bio- 
mathematics and in social and engineering sciences - require 
the clarification of at least the nature of the behaviour of 
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solutions of these equations in the entire domain of their 
existence without explicitly solving them. In this direction, 
the development of the qualitative theory owes its legacy 
to Poincare and Liapunov towards the end of last century/ 
forms a significant contribution. The basic problems of 
study of this theory are (i) the existential analysis and 
( ii) the properties of solutions such as positivity/ bounded- 
ness and stability. 

The qualitative theory of integral equations largely 
owes its origin to Vito Volterra and Ivar Fredholm. In 1896 
Vol terra encountered integro-differential equations while 
investigating the dynamical processes mainly in biology. 

With the advent of Fredholm theory in 1900 and with its 
tremendous Impact on mathematical and functional/ analyst?S/ 
the study of Volteirra integral equations was relegated into 
the background for quite some time* However/ during the 
last thirty years Volterra equations have emerged vigourously 
in many applied fields such as nuclear reactor dynamics (cf. 
[28])/ autcxnatic control systans (cf, [s] /[9] /[l2] / [26] ) / 
the theory of industrial inventory (cf . [38 ] ) , and in biology 
and hereditary mechanics. 

Classically there are three methods available for 
the qualitative study of these equations viz,/ Perron's 
technique, Liapunov direct method and Popxjv's method. 

Perron's technique basically involves the study of perturbed 
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linear systems -v^en the perturbation functions are sufficiently 
small in a certain sense. This technicfue has been onployed 
extensively by various authors (cf, [16] , [17] / [33] / [35] / £36] ^ 
[ 37 ] ) for Vol terra integro-differential equations via the 
resolvent function and the variation of constants formula. 

In Popov's method/ the integral transforms are used and as 
such it is severely restricted to integral equations of 
convolution type, A common feature of the systems to viiich 
Popov's method could be applied is that the kernel is 
integrable. This method has been exploited by various 
authors (cf. [s], [ 9 ] , [I 2 ] / tl3] / [si] ) to study the stability- 
properties of solutions of Volterra integral equations, 
Liapunov direct method in-volves the construction of a scalar 
function, -the time derivative of which should satisfy a 
negative estiraa-tse along the solutions of "the equa-tion. For 
integro-differential equations, Liapunov functionals are 
constructed instead of functions and perhaps this is na-fcural 
to expact due to the hereditary nature of Volterra equations. 

A number of results (cf. [1] jr[2 ], [25]) have been reported and 
a general result due to Levin involving non— con-volu-tion 
kernels could be found in [13, chapter 5.1], In recent 
years Liapunov— Raxumikhin technique (cf, [41], £21]) vjhich 
involves functions rather than func-tionals has been used -to 
study integro-differential equations by many authors (cf, [15] , 
[24], [ 47 ], [48] ) and a good number of results have been 
accumulated. In addition -tx* -the above techniques, -the ■t±ieory 
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of ac3itiissibilit.y developed by Massera and Schaffer for 
ordinary linear differential eqaations has been further 
extended in recent years for Vol terra integral ecjuations 
(cf. [I0]^[ll]^[l3] , [ 34 ] ,[ 40 ]). 

1 .2 Brief Review 

The qualitative properties of solutions of Volterra 
integro-dif ferential equations have been the focus of study 
in many works [cf. 1-3,12/15-18,24-26,28,29,31,35,37,42,43/ 
47—50] under different hypotheses and various results have 
been acciamulated. 

Friedman [l4], Kornienko [ 22 I , Levin [27], Ling [30], 
Miller [33] and Padmavalli [38] among others have studied 
the positivity of solutions of Volterra integral equations. 

While non-convolution equations are considered in [22] , the 
convolution equations are discussed by several others [14,27, 
30,33,38] » For integral equations of the form 

t 

x(t) sr f(t) - I a( t,s) x(s)ds, O < t < “> (1.2.1) 

o 

an application of Neumann series yields the positivity of 
solutions of ( 1 . 2 . 1 ) if the kernel a(t, s) is negative and 
the source function f(t) is positive. Ling [30] has 
discussed the positivity of solutions of (1.2.1) in the 
convolution case x^en the kernel a(t) is positive and its 
derivative a'^(t) is negative. However these results are not 
in general applicable to Volterra integro-dif ferential equations. 
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Therefore it is interesting to investigate positivity of 
solutions of 

t 

x*{t.) = f(t)-/ a( t, s) x(s)ds/ x(0) = ® ^ (1«2*2) 

o ” 

l^articularly when the kernel a(t^s) is positive. An attempt 
has been made in this direction for linear and non-linear 
Volterra integro-differential equations. 

When zero is not a solution of the system of ordinary 

differential equations, eventual asymptotic stability of the 

set £0} is discussed in [23,51]. Very often the source 

function f(t) is not zero in (1.2.2) and it is therefore 

appropriate to study the eventual properties of the set £0} 

for (1.2.2). As such the eventual uniform asymptotic stability 

of the set £0} for linear and non-linear Volterra integro- 

differential equations are discussed here using resolvent 

functions and variation of constants formula* Quite recently 

Seifert [47,48], Grimmer and Seifert [15], and Lakshmikantham 

and Rcmia Mohana Rao [24] have studied stability, boundedness, 

uniform stability and uniform asymptotic stability of Volterra 

integro-differential equations by employing Liapunov- 

Raiumikhin technique which involves Liapunov functicffis rather 

than functionals. Hiis technique is e:xploited further in this 

2 

thesis to discuss L —stability of Volterra integro-differential 
systons* 

The study of linear Volterra integro-differential 
systems (VIDS) very often rec[uires that the coefficient matrix 
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be stable. However the notable exceptions that have dispensed 
with this condition on the coefficient matrix have been the 
works of Levin [25] ^ Grossman and Miller [17] » Grimmer 
and Seifert [15] and Burton [3]» By adopting an entirely 
different approach involving an arbitrary function the 
exponential stability properties of solutions of Volterra 
integro-diff erential eguations# vhen the coefficient matrix 
is not necessarily stable, are investigated in the later part 
of the thesis, 

1.3 Layout of the Ihesis 

The present thesis attempts to study some of the 

qualitative properties of solutions of Volterra integro— 

differential equations such as positivity, boundedness, 

eventual uniform stability, eventual uniform asymptotic 
2 

stability and L -stability. Dichotomy properties of these 
equations are also discussed. Ihe entire work of this thesis 
is divided into five chapters. 

The first chapter deals with historical survey and a 
brief review of the work. Various techniques available in 
the literature and recent contributions towards the gradual 
development of the theory are indicated. 

Chapter 2 deals with the positivity and boundedness 
of solutions of Volterra integro-differential equations. 

It is clear that one can assert from Neumann series that the 
solutions of (1,2,2) are positive if the kernel 'a' is 



7 


negative and the weight function * t* is positive. However 
the same cannot be concluded if 'a' is positive. An attempt 
has been in this chapter to ascertain the positivity of the 
solutions of (1.2.2) particularly when 'a^ is positive. 

The technique employed in this chapter is to transform (1,2.2) 
into an equivalent equation whose kernel would be negative by 
a proper choice of the auxiliary function. The results 
obtained here are compared with those of earlier authors and 
examples are provided to illustrate the effectiveness of the 
approach. This study includes both convolution and non- 
convolution equations. The case where the solutions exist 
over a finite time interval is also discussed. 

Motivated by the fact that zero is not a solution of 
(1.2,2) when the forcing function is different from zeix>/ 
eventual uniform/ stability and asymptotic stability properties 
of solutions of (1.2 .2) are investigated in Chapter 3, The 
main tools in this chapter are the resolvent kernels and 
variation of constants formula. Perturbed equations are also 
studied under a more general class of non-linear perturbations. 

In Chapter 4, Liapunov-Razumikhin technique is employed 
2 

to Study L -stability properties of solutions of Vol terra 
integro-differential systems. In order to proceed with the 
avowed intention/ a more general result on L^- stability of 
solutions of a functional differential systan of delay type 
vdiere the delay becomes unbounded as t -► <»/ is derived and 
thoi applied to Vol terra integro-differential systems. 
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Finally In Chapter 5 the exponential asymptotic 
stability properties of solutions of linear Vol terra integro- 
differential systsns where the coefficient matrix is not 
necessarily stable/ are obtained. Dichotomy properties of 
these systems also find a place in this chapter. 



CHAPTER - 2 


POSITIVITY AND BOUNDEDNESS OF SOLUTIONS 
2.1 Introduction 

Since long there has been a persistent effort with 
quite an amount of success to descem various properties - 
like boundedness / stability, asymptotic behaviour, positiveness. 
Square integrability and so on - of the solutions of integral 
equations of Volterra type that occur in a variety of problems 
which are heriditary in nature such as in nuclear reactor 
dynamics, biological models etc. In this chapter, positivity 
and boundedness of solutions of Volterra Integro-differential 
equations are studied, 

Friedman [14] , Kcsnlenko [22]/ Levin [27], Rina 
Ling [30], Miller [33], Padmavalli [38] and Rama Mohana Rao 
and Srinivas [45] among others have studied the positivity 
and boundedness of solutions of Volterra integral equations 
and many interesting results have been reported. While 
nonconvolution kernels are taken into account by one of 
earlier workers [22], the kernels considered ty several 
others [14,27,30,38] are of convolution type. The results 
obtained in [45] constitute the main core of this chapter. 

In the case of integral equations of the form 

t 

x(t) = f(t) - / a(t,s)x{s)ds/ 0 < t < «» (2,1,1) 

o "* 



one can assert by virtue of Neumann series (cf. [52]) that 
if the kernel aCt^s) is negative and the source function f(t) 

IS positive/ then the solutions of (2*1«1) are positive in 
the interval of their existence. It is interesting to probe 
the conditions on fCt) in order to obtain the positivity of 
solutions of integro-differential equation 

t 

x'(t) =s f(t) - / a(t/S)x(s)dS/ x(0) = X > 0 (2*1 *2) 

o ° 

where 0 < t < «>/ especially vdien the kernel a(t/s) is positive* 
The objective of this chapter is to investigate sufficient 
conditions for positivity and boundedness of solutions of 
Vol terra linear and nonlinear integro-differential equations 
of both convolution and nonconvolution type* Our approach 
herein involves the construction of an auxiliary function 
which/ in certain cases, turns out to be the usual resolvent 
kernel* The advantage of such an approach is to transfer the 
( i) responsibility of the original kernel to a better kernel 
and (ii) unwanted behaviour (if any) of the original kernel 
to the source function, through a proper choice of an 
auxiliary function* Examples and remarks are provided at 
appropriate places to illustrate the generality of the 
conditions and the fruitfulness of the results. 

Under suitable conditions on a(t, s) (for exait^le, a(t,s) is 
absolutely continuous on every ccai^act subset of HR x 3R) , the 
equation (2*1,2) can be written as 
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t t t 

x'Ct) = f(t.)- 3 ?^ / a(t;,s)ds-/ [/ a(t,u)da] x''Cs')ds,(2-1.3) 

o os 

t 

Therefore, if f(t) - / a(t, s)ds is positive and 

t 

/ a(t, s)ds is negative, then in view of Neumann series, one 
s 

can conclude that x'(t) in (2.1 *3) is positive and hence the 
solution x(t) of (2.1*2) is positive as > 0 (such a 
conclusion is quite obvious if a(t,s) is negative and f(t) 
is positive) • This fact has been extensively used in our 
subseguent analysis. 

While Section 2.2 is devoted to basic results, 
sections 2.3 and 2.4 cover the linear non-convolution and 
convolution equations respectively. Non-linear equations 
are discussed in Section 2.5. Finally, in Section 2.6 
similar results over a finite interval are obtained which 
assert that an heriditary system won't go 'critical' in that 
interval * 

2.2 Basic Results 

In this section we shall derive some basic results 
which will provide the necessary leverage to workout the 
planned trade off between the kernel a(t, s) and the source 
function f(t) . 

Lanma 2.2.1 If a(t, s) is continuous onO<s<t<<o and 
f(t) is continuous on 0 < t < «», then the eqpaation (2.1.2) is 
equivalent to 
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t 

y'(t) = h(t)-/ b(t,s)y(s)ds^y(0) =XQ>0,0<t<«> (2.2*1) 

o 

where 

b(t, s) = a(t, s)+cpg(t/s)-/ 9(tyti)a(u/s)du 

® ^ ( 2 , 2 . 2 ) 

+ g>(t,t)/ a(u,s)du, 
s 

<j>(tys) being a c^-function onO<s<t<'» 
and 

t 

h(t) = f(t)+XQ <?( t,t)- 3 ^ 9 ( t,0) +q)( t,t) / f(s)ds 

° (2.2.3) 

t 

- / ^(t,s) f(s)ds. 

o 

Proof. Let x(t) be a solution of (2.1.2) existing on the 
interval O < t < «>. Consider the identity 
t 

/ ^-( t,s)x(s)ds = <p(t/t)x(t) - ^(t,0) x^ 

° t- 

- / <p(t,s) x'(s)ds. 

o 

Substituting x(t) and x'^(t) from (2.1.2) and using Fubinis 
theoroTi/ we get 
t 

/ (Pg( t/s) x(s)ds = x^ <P(t^t) - Q)(t,0) 

o 

t t ' 

+ <(p(t,t)/ f(s)ds-/ ^(t,s)f(s)ds 
o o 

t t 

- (})(t/t) / [/ a(s^T)ds] xCt) dT 

o T 

t t 

+ / [/ a(s,T) 9 {t^s) ds] xC'ddT (2.2*4) 

o T 
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Therefore, it follows frcxn (2.1*2) , (2.2.2), (2.2.3) and (2.2.4) 
that 

t 

/ b(t,s) x(s)ds = ~x'(t) + h(t), x(0) = Xq 
o 

This shows that every solution of (2.1.2) is a solution of 

(2.2.1) . For converse, let Y(t) be a solution of (2.2.1). 

Then by taking <ji(t, s) = 1, we have from (2.2.2) and (2.2.3), 
b(t,s) = a(t,s) and h(t) = f(t) and hence y(t) satisfies 

( 2 . 1 . 2 ) . 

Remark 2.2.1. If we choose b(t,s) = O i.e. 9 (t,s) is a 
solution of (2.2.2) with b{t,s) = 0, then frcm (2.2.1), the 
solution x(t) of (2.1.2) can be obtained on integrating (2.2.3) 
between 0 and t. 

A result similar to Lemma 2.2.1 for convolution kernels 
is given below. 

Lsnma 2.2 .2. If a(t) and f(t) are continuous on O < t < «>, 
then the following equations are equivalent : 

t 

x'(t) = f(t) - / a( t-s)x(s)ds, x(0) = (2.2.5) 

o 

t 

y'(t) = h(t) - / b( t-*s) y( s) ds , y(0) = x (2.2.6) 

o ° 

where O < t < «>, and 

t" t 

b(t) = a( t) +<p ' ( t) +/ q5(t->T)a(T)dc-<p(0) / a(T)dt, {2.2.7) 

o o 

<p e c^[o,oo)^ 
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t 

h(t> = X <p(t)-x (?(0) + / <p(t-T)f(T)dT 

O 

t 

- <p(0) S f('c)dT: + f(t) (2.2.8) 

o 

The following result is somewhat close to the usual 
resolvent equations (cf. [16] and [37} ) and it is useful in 
the study of nonlinear Volterra integro-diff erential equations 
(see section 5). 

Lemma 2.2.3. If a(t) and f(t) are continuous on 0 < t < 
then every solution x(t) of (2.2.5) satisfies 

t 

x(t) =f(t) — / c(t-T)x(T) dT (2.2.9) 

o 

on O < t < <», v^ere 

t 

c(t) = (p'(t) + / q5(t-T)a(T)dT , (2.2.10) 

o 

t 

f(t) = x^ «p(t) + / f(s) <P(t-s) ds^ (2.2.11) 

o 

and 9 e c [o,<»),?(0) =1. 

Remark 2.2.2. Choose <|)(t) (for example, ^(t) = e + te 
with a(t) = such that c(t) = 0. Then (2.2.10) takes 

the form 

t * 

<p'(t) = -/ <P (t-T)a(T)dT ^ Q5(0) = 1 
o 

Thus 9 (t) would agree with the differential resolvent to the 
kernel a(t) and the eq[uation (2.2.9) gives the solution x(t) 
of (2.1.2) on O < t < «> in terms of the resolvent (pit) (see 
Grossman and Miller [1^ , equations (A) and (4) with A(t) = 0) . 
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Remark 2.2.3. 


It should be noted that the above lemmas are true for 
solutions defined over a finite inteirs/al J = [O/Ot] with 
appropriate changes. 

2.3 Non Convolution Equations 

In this section the positivity and boundedness 
properties of solutions of (2.1.2) are studied through a 
proper choice of the auxiliary function ^(t, s). 

Theorem 2.3.1, Suppose for 0 < s < t < (i) a(t/s) >0 

t 

( ii) X(t) = f(t) - y:^ / a(t^u)du is an increasing function 

t t 

of t and (iii) / / a(t,T) dx ds < 1 hold. If xCt) is any 

o s 

solution of (2.1.2) on , then the following estimates 

are satisfied on [o^«>) : 

(a) 3^ < x(t) < 3^ + P(t) 

(b) 0 < x'(t) < f(t) 

t 

where F(t) = / f(s)ds, 

o- ^ 

Proof* Take <()(t,s) = / a(t#u)du. Then <j)(t,t) = 0, 

s 

Since a(tys) > 0/ it follows from (2.2.2) that 


t t 

b(t, s) sa a(t/ s) -a(t, s)-/ a(u^s)(/ a(t,x)<3x)du 

s u 

< 0 . 

Further, using (2. 2.2) and (2.2.3) we obtain 



16 


t t t t 

b{t/u)cfc. = f(t)+x^/ [/ a(u/s)(/ a(t^T)dT)<3uj ds 
o o s u 

t t "t 

-/ (/ a( t/u) du) f( s) ds-x^ S a(t/u)du 
os o 


t 

= f(t) - Xq / a(t,u) da 
o 

t u t 

+ x^J is a(u<s)ds)(/ a(t,T)dT)du 
°o o u 

t t 

- / (/ a(t,u)da) f(s)ds 
o s 

t 

= fCt) - Xq / aCt^u) da 
o 

t t s 

-/ (/ a(t,x)dT) [f(s)-x^/ a(s,u)du] ds 
os o 


t: t 

= Kit) ~ S [/ a(t,T)d'r] X(s) ds 
o s 

t 

Thus the assumptions (ii) and (iii) yield hCt)-x^/ b(t, s)ds > 0< 

o 

Hence, in view of the observations made toward the end of 
section 1, it is clear that the solution y(t) of (2*2«1) 
together with its derivative, is positive for all t e [o^<») . 
Moreover, we have O < y^(t) < h(t) on O ^ t < «»• CBierefore^ 
by invoking Lemma 2«2»1, we conclude that the inecjualities 
(a) and (b) hold for all t e [ 0 ,oo) and this completes the 
proof. 


Corollary 2.3.1. If, in addition to the assumptions (i), 

1 

(ii) and (iii) of Theorem 2.3.1, suppose f e L [0,oo) , then 
the solution x(t) of (2.1.2) is positive and bounded on [0,«.) . 
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Example 2.3,1, Choose a{-t,s) = e /(1+s) for O £ s < t < «»• 
It can be seen that the conditions ( i) and { iii) of TheorCTi 
2.3,1 are satisfied. If we take f(t) = 3-e~^ and x^ = 1/ then 
the condition ( ii) also holds. 


Remark 2,3,1, Let 


and 


t 

Pj^(t) = f(t) - Xq / aCt/u)du 

o 


t 

KCt^s) = — / a(t,u)du 
s 

for 0 < s < t < <», Then the equation (2.1.2) reduces, upon 
integrating by parts (see equation (2.1,3), to 


t 

x'(t) = P 3 _(t) + / K(t,s) x'(s) ds (2.3.1) 

o 

on 0 < s < t < «», In order to obtain an estimate similar 
to (b) of Theoran 2,3,1 on the solution x'(t) of (2,3,1), 
Kornienko [22 , Corollary 2 and remark 4] assumed the 
condition 


If we choose a(t, s) = e /(1+s) and f(t) = 3-e for 
0 < s < t < “>, and = 1, thoi 


[§r ( | -' r£ f^)] - - 0.4975 

Ldt P^^t; -"^=0.001 “ 

SssO m OOOl 

Thus the condition 3,2 of Kornienko [22] does not hold. 
However, all the assumptions of Theorem 2,3,1 are satisfied 
(see Exan^le 2,3,1) • 
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2»4 Convolution Equations 


We shall now give a result similar to Theorem 2. 3*1 
for convolution equation (2.2.5). 


Theorem 2.4.1. Suppose for O < t < «> ( i) a(t) > O 
(ii) a'(t) < 0 Ciii) f(t) > O, f'(0) > a(0), f'^Ct) > a(0)f(t) 

hold. Then if x(t) is a solution of (2.2.5) on 0 < t < «», 
then we have. 


(a) < x(t) < Xq + F(t) 

(b) O < x'(t) < f(t) 

t 

where P(t) = / f(s)ds. 


Proof. Take <p(t) = — / a(T)dx . Ihen from assumptions 

o 

( i) and (ii) it is clear that Q)(0) =0, <P^(t) = -a(t) < O 
and (p"(t) = -a'(t) > O. Moreover, from (2.2.7) and (2.2.8) 
we have 

t t--T 

b(t) = -/ [ / a(s)ds] a(x)dT < O for O < t < «>, 
o o 

h(0) = f(0) > O and 

t 

h'(t) = XQ 9 '(t)+/ g)'(t-x)f(T)d'r + f'(t) (2.4.1) 

o 

Using the assumption (iii) and (2.4,1), we obtain 


h'(0) = Xq (P'CO) + f'(0) = f'(0) - x^ a(0) > 0 

and 

t" 

h'^t) = Xq t) +< p' ( O) f ( t) +/ t-x) f(T)dx + f^^t) 

o 


> f"(t) - a(0) f(t) > O. 
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Thus we have h(t) >0 and bCt) < O for all t 6 • Hence 

the Lenma 2 *2 *2 together with the obsearvation at the end of 
o action yield that tbe solution x(t) of (2«2*5) along with 

its derivative x*(t) are positive on • Since x(t) >0 

and a(t) > 0, it is clear from (2»2»5) that 0 < x^ ( t) < f(t5 
and hence (a) and (b) are satisfied for all t 6 • 

Remark 2.4.1. If, in addition to the asstimptions of 
Theorem 2.4.1, suppose f(t) e L , then the solution 

x(t) of (2-2.5) is also bounded on [0,<»). 

Remark 2.4.2. While discussing similar properties for 

solutions of Vol terra integral equations of convolution type. 

Ling [30], Miller [33], and Levin [27] have assumed that 

a(t) >0 and a'(t) < O. When we integrate(2 .2 .5) between 

0 and t, then the resulting integral equation has the kernel 
t 

/ a(s)ds and hence the above assumptions cannot be satisfied 
for such a kernel. Therefore the results of [27], [30] and 
[33] are not applicable for (2.2 oS). Our endeavour here has 
been to find conditions such that the auxiliary (or resolvent) 
kernel b(t) is negative so that the solutions could be positive. 
More specially, our results assert that the resolvent could 
be negative without the kernel being negative. 

2.5 Nonlinear Equations 

Having known the conditions on the kernel a(t) and 
the source function f(t), for the solution x(t) of the linear 
convolution equation (2.2.5), we shall now discuss similar 
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properties of solutions of the nonlinear equation 

t 

x'(t) = fCt)-/ aC t-s)x( s)ds+G(t/x(t) ) #x(0) = x > O (2*5. 1) 

o 

o 

where 0 < t < oo under growth conditions on the nonlinear 
perturbation function G(t,(f)) . Assume that 

(Hj_) : GCt,g>) > O for all t > O and «p > 0 
(H 2 ) ; there exists a positive continuous function 
X(t) such that 

lG(t,<p) 1 < X(t) (1+1 cpi) 

where lim X(t) = 0/ 0 ^ t < <»« 

(H^) : the resolvent <?p(t) associated with the kernel 

a(t) is positive for all t e [o^oo) and 
1 

g)(t) e L [o,<») . 

Theorem 2,5«1. Suppose the hypothesis ^ and (H^) 

are satisfied* If the solution yCt) of (2 .2 *5) with y(0) = x^ 
is positive and bounded on [o^«>), then every solution x(t) 
of (2.5.1) with x(0) = Xq satisfies the inequality 

Xj^ < x(t) < 
for all 0 < t < »• 

Proof. From Lemma (2.2.3), Remark (2.2.2) and (2.5 J.), it 
follows that 

t 

x(t) =s y(t) + / <p(t-s)G(s/x(s) ) ds. 

o 
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Define 

H(s,x) = G(s,x^) for s e [o^<»)/ x b 3R (2.5.2, 

where = x if x > O 

= 0 if X < O. 

Thus H(s,x) > O for s e [o,co). x eB, Let p(t) be 
a solution of 

t ” ^ 

p(t) = y(t) + / q)(t-s) H(s,pCs)) ds (2.5.3) 

o 

Therefore# if y(t) is positive on [o#<») , then (Hj^) # ^^ 3 ^' 
( 2 . 5 . 2 ) and ( 2 . 5 . 3 ) imply that p(t) is positive on [o,oo) 
and hence the solution x(t) of (2.5.1) is positive on [o#<») . 
Moreover# from the boundedness of y( t) on [0/'») # the 
CoTOllary 2.1 of Nohel [37] yields that x(t) is also 
bounded on [o,oo) . This completes the proof. 

2.6 Behaviour of Solutions over a Finite Interval 

This section deals with the positivity of solutions 
defined over a finite interval# of Volterra integro- 
differential equations of both convolution and nonconvolution 
type. The generality of the conditions and the fruitfulness 
of the results are illustrated by means of examples. 

Theorem 2.6.1 
Suppose 

(i) O < a(t) < cf# 0 being a constant 

and 


(ii) f(t) >0 hold. 
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Then there exists an interval [o^aj / where a is such that 
2 

(aa +ba)<2/bisa positive real number/ on which every 
solution of (2.2.5) satisfies 


0 < x(t) i Xq + F(t) 
t 

where P(t) = / f(s)ds. 
o 

Proof. From (2.2.11) we have 

t 

f(t) = Xq c?(t) + / f(s) (p(t-s)ds. 

o 

2 

For a given ct > O there exists a b > O such that 0(X + bot < 2. 

Choose <P(t) = - ^t^ — ^ t+1 for all t e [O/tt) , 

Then <p(t) >0 [0,a] and (p(0) =1. Also <p'( t)=^ {t)-b/2 ^ C 
on t s [O/d) • 

Thus from (ii) we have 5f(t) > O on [ O/dj • 

Equation (2.2.10) gives 

t 

c(t) = (p'(o) + / (p(t-T:)a(i:)dT 

o 

This implies that c(0) = 9^(0) = -b/2 < 0. 


Therefore from (i) and the fact q)'(o) < 0 for all t > 0/ 
it follows that 


t 

c' (t) = ^"(t) + cq)(0)a(t) + / q)'(t-T)a(x)dT 


< -a + a + / <p'(t-ir)a(T)dT 


o 


u:. 


< O 
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Thus we obtained c(o5 < 0 / c' Ct) < O and hence c(t) < O, 

Now fzxjm the observations made towards the end of section 2«1^ 
WCt) > 0/ and c(t) < 0 , every solution of ( 2 # 2 .9) is positive. 
The application of Larona 2.2.3 now yields that every solution 
x{t) of (2.2.5) is positive. 

Since x(t) is positive on [o,a] and a(t) > 0 for 
t e "[o^a] we have from equation (2.2.5) that 

x''( t) < f(t) . 

This in turn implies that 0 < x(t) < x^ + F(t) for all 
t e [O/tt] , 

Theorem 2.6.2 

Assxime that 

(hj^) : for 0 < s < t < T 

(i) a(t,s) > O, < a(s,s) < 

(ii) a^(t,s) < 0 

(iii) O < ag(t,s) < M 3 

(h 2 ) : for a given a(0 < a < T) there exists positive 
constants 6 and K such that 

(i) p > ra + 6 a^ 

(ii) 0 M 2 

(iii) M 3 < XM^ 


and 
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(h3) : 

(i) f(0} > 0 

(ii) f'(t) + 5 . g? .) f(t) 

(5a +m-3) 

Xq 3 x(r+ 26a) 

> — 

(0-ya-6a ) 

a,3,ry6/ and X being the same constants 
as in (h2) • 

Then on the interval [o a] every solution x(t) of ( 2 » 1 « 2 ) 
satisfies the following inequalities 

(I) O < x(t) < 3 ^ + F(t) 

and 

(II) 0 < x'(t) < f(t) 

t 

where > 0 and F(t) = / f(s)ds. 

o 

Proof* For O < s < t < a (fixed a < t), define 


<}>(t/s) = p(t) q(s) 


where p(t) s s- and q(t) = (g - yt - 6‘t ) • 

-p+yt +6t 

By (h 2 ) (i), q(t) > O, q'(t) < O and q*(t) < O and t S [0,a] 

( 2 . 6 . 1 ) 


Also p(t) > O and p'(t) 


X(yt26t) 

(-p+ya+6<i^)^ 


> O for t c Co,a] 

(2.6.2) 


Moreover ^P(tj,t) = p(t)q(t) =‘K > 0 
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From (2.2.2) and (h 2 )(ii) we have 

b(0*0) = a(0,0) + pCO) q'(0) < - ^ < O 

Differentiatxng (2.2.2) with respect to ^t' we get 

t 

b^(t,s) = aL{.( t^s)+p'( t) q' (s)-/ p'(t)q(u>a(u/s)du. 

s 

From (hj^) ( i) , ( ii) (2.6.1) and (2.6,2) it follows that 

b^(t,s) < O on t e [o,a] . 

Differentiating (2.2,2) v;ith respect to 's^ we get 

t 

bg(t^s) = ag(t,s)+p(t)q"'(s)+p(t)q(s)a(s/s)-/ p(t) q(u) ag(u ,s)cu 

s 

t 

- \a(s,s) + X / ag(u^s)du (2.6,3) 

s 

Since p(t) >0, q^(b) < O and ag(txs) > O we have 
t t 

X S a (u,s)du-/ p(t)q(u)a (u, s)du < O for all 0 s < t < a 
s s 

(2.6.4) 

Hypothesis (hj^) / (h 2 )(iii) now gives 

ag(t/.s) -Xa(sj.s) < "*XH 2 < O (2.6,5) 

From (1^) (ii) we have 

O > 0M2-26 > 3a(s,s)-26 > (0--rs-6s^)a(s,s)-2 6 (2.6.6) 

Further (2.6.6) together with p(t) >0 gives 


p(t)q^(s) + p(t)q(s)a(s/s) < O 


(2.6.7) 
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Now substituting (2.6.4)^ (2.6.5), (2.6.2) in (2.6.3) we get 

bg(t, s) < O for O S s 5. t < a. 

Since b(0,0) < 0, b^(t,s) < 0 and bg(t,s) < O, it is clear that 
b(t,s) < 0. 

Prom (2.2.3), and the fact p'(t) > 0, we have 

h(0) = f(0) > 0 

and 

t t 

h(t) > f (t) + xx^-pXQp(t) + X / f(s)ds-p(a)/ q(s)f(s)ds 

o o 

( 2 . 6 . 8 ) 

Define the right side of the inequality (2.6.8) by x(t), 
then X(0) > O and 


x' ( t) = f* i t) -pXj^p' ( t) + X f ( t) - 


Moreover from (h^) ( ii) we have 


3-rcx-6a' 


^q(t) f (t) 


(2.6.9) 


f'(t) +Xf(t) + 


J2l 


f(t) > 


x^8x(r+26a) XQ^^(y+26t) 

27T ^ 77~im2T2 


(6a^+<5a-0) (8-ya-6a ) (0-rt-6t^) 


This implies that 


f'(t)+Xf(t) > PXq p'(t) 1 


§2l 


— f(t) > p'(t) 

(6a^+ra-3) 

(6a^+y®-*^) 


> 33{^ p'(t) + 


(3-ya-5a ) 


q(t)f(t). (2.6.10) 
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Therefore (2.6,9) and (2.6.10) give X^(t) > O for t e [O/a], 
This together with X(0) > 0 shows that X(t) > O for all 
t S' [C,cc,] , Ihus from (2.6.8) we have h(t) > X(t) > O. 

Finally the observe non s made towards the end of section 
of 2.1 ensures that all the solutions of (2.2.1) together with 
their derivatives are positive. Therefore in view of Lannma 
2,2.1/ the assertion of the theorem follows. 

Ranak 2.6,1 

In theorem 2.6,2 the hypothesis (h^) on f can be 
replaced by 

ih*) : (i) f(0) > X( 2 ^ - 1) 

j3-ra-6a 

(ii) f^(t) > (I5L±.^) f(t). 

Example 2.6,1 

Select a(t,s) = for 0< s< t<a (=1) and 35^=1 

Then a(t/s) > O and a(S/s) =1, Therefore = ^2 “ 

Moreover a^(t/s) = — e < 0 and ag(t/s) = e > 0 and 

ag(t/s) < 1 (= M^) . 

Let ^ = 3/2/ y = 1/3/ 6 = 1 and X = 6 so that (h^) of 
theorem 2.6,2 is satisfied. 

Choose f(t) = then f(o) = 1 and 
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f'(t) 4.M21|±S2dL f(t) > I000e^°°°^ - > 952 

(sa'^+ra-g) 

X ex (y+26a) 

> 756 = -2 . 

(3-ya-6a ) 

Thus all the conditions of theorem 2.6.2 are verified. 


Example 2.6.2 


If we taka f(t) = 50e5°t i„ example 2 .6. 1, then 
(h|) holds. However, (h, ) does not hold. Therefore, the 
oondrtions(h 3 } and (h|) are some what independent. 



CHAPTER - 3 


EVENTUAL STABILITY PROPERTIES OF SOLUTIONS 
3.1 Introduction 

The purpose of this chapter is to study the eventual 
uniform stability and eventual uniform asymptotic stability 
of set {0} for the integro-differential equation 

t 

x'(t) = f(t) + / a(t,s)x(s) ds (3.1.1) 

o 

and the peirturbed equation 

t 

z'(t) = f(t)+I a( t, s){ z(s) +g(s *’2 Cs)JL}-ds. (3.1.2) 
o 

Eventual uniform asymptotic stability is a natural generali- 
sation of uniform asymptotic stability in which it is not 
assTimed that the zero function is a solution of (3.1.1), 

Infact, this would happen if f(t) ^ 0 in (3.1.1). The main 
tools used are the resolvent function and the variation of 
constants formula. 

In the first part of this chapter, we shall investigate 
sufficient conditions for eventual uniform asyn^totic stability 
of the set {0} for (3.1.1) under certain hypotheses on the 
resolvent function r(t, s) (corresponding to a{t/s)) and the 
weight function f(t). Convolution equations are also discussec... 
In many practical problems if {0} is eventually uniformly 
asynqptotically stable for (3.1.1) then it is interesting to 
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know Idle class of pertixrbation functions g which preserve 
the same property for (3.1.2). In the later part of the 
chapter an attempt has been made in this direction under 
certain growth conditions on the nonlinear perturbation 
function g, 

3.2 Definitions and Basic Results 

The following definitions are useful in our subsequent 
discussion. Let zCt/t^/zCt^)) be any solution of (3.1.2) with 
z( t^y t^ 2 (tQ) ) = 2 (tQ) existing for all t > t^. 

Definition 3.2.1. 

The set {0} is eventually uniformly stable (EvUS) if 
for every e > 0, there exists a = a(e) > O and 6=6(s) >0 
such that 

1 2 (t/tQ# 2 (tQ) ) I < s for all |z(t^)l < 6 and t > t^ > a. 
Definition 3.2.2. 

The set {0} is eventually uniformly attracting (EvUA) 
if there exists 6^ > 0 and > O and if for every h > O 
there exists T = T(T)) > O such that 

I z( t/tQ,z(tQ) ) I <'V for lz(tQ)l < d^^t^ > and t > t^+T 
Definition 3.2.3. 

The set £0} is eventually uniformly asymptotically 
stable (EvUAS) if it is both E\^S and EvUA. 
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Definition 3.2*4. 

Let be the class of monotonic, non— negative 
functions X ( . ) defined on [ 0,«>) such that X(P) - Oas p-*0. 

Definition 3.2.5. 

Let be the class of monotonic, non— negative 
functions v(,) defined on [0^<») such that v(t) -* 0 as t -• <». 

The definitions 3.2.1 - 3.2.3 are equally applicable to 
the convolution equation 

t 

y'(t) = f(t) + / b(t-s)y(s)ds (3.2.1) 

o 

and also to the nonconvolution equation (3.1.1). 

We need the following basic result due to Grossmaui and 
Miller [16] expressing the solution x(t) of (3.1.1) in terms 
of the resolvent function r(t/s) satisfying the equation 

t 

~ r(t/s) = -/ r(t^u) a(u/s)duy r(t/t) = 1 (3.2.2) 
s 

on the interval O £ s £ t. 

Let LL^(S) be the set of all measurable functions on a 
set S such that the seminorms are finite for all ccxnpact subsets 
of S, 

Lemma 3.2.1 ([l6]) 

Assume a(t,s) e LL‘( H x »*’) . Then the function r(t/s) 
satisfying the equation (3.2.2) exists on O < s < t and is 

continuous in (t^s). exists a.e. on O < s < t^ is in 
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LL^( 3R^ , Moreover^ given any and a function 
f(t) e LL^( 3R*') equation (3.1.1) iS equivalent to 

__ t 

x{t) = r{t/0) ^ + / rCt^s) f(s)ds. (3.2.3) 

o 

For the convolution equation (3,2.1) if c(t) is the 
resolvent function satisfying the equation 

t 

c'(t) = / b( t-s)c( s)ds/ c(0) = 1 (3.2.4) 

o 

then the solution y(t) of (3.2.1) with y(0) = y^ is given by 

__ t 

y(t) = c(t) y_ + / c(t-s) f(s)ds. (3.2.5) 

° o 

We shall now state a result for vol terra integro- 
differential equations which is analogous to lerama 3.1 of 
Strauss and Yorke [51]* 

Lemma 3.2.2* 

Let £0} be EvUAS for (3.1.1). Suppose there exists 
X e and ^ and for some r^ > 0/ each ra e C0,rj^), 

and each x > o, there exists F_ _ e K such that for each 

T fin ^ 

e (O/T] , each t^ > a(rj^) and each solution z(.) of (3.1.2) 
s atis f ying 

lz(tQ)l < 6(rj^) and m < lz(t)| < r^^ for t^ < t < t^+ t j 
there exists a solution x(t) of (3.1.1) such that 


|z(t)-x(t)l < X( 1 z(tQ) I ) 3^,(-t:-t;o) + 
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for all < t < Ihen {0} is EvUAS for (3,1.2). 

To proof IS virtually similar to that of Lemma 3*1 of 
Strauss and Yorke [Sl] and hence onitted, 

3,3 Eventual Uniform Stability 

In this section we shall investigate sufficient 
conditions for the eventual uniform stability of the Vol terra 
lntegro~dlff erentlal equation (3,1,1), A few propositi a us have 
been proved first which assert the Integrability and boundedness 
properties of solutions of (3.1.1), These results have been 
used to obtain eventual uniform stability of the set {O}. 

We shall now give an account of basic hypotheses on 
the resolvent function r(t,s) and the weight function f(t). 

The generality of the hypotheses and the fruitfulness of the 
results are Indicated through various remarks and examples. 


Assume that 

(H^) t Ir(t^s) 1 < 

(H 2 ) : f(t) e [ 0 ,co) 


where is a positive constant^ 

CO 

(l.e. / lf(s)lds < M^/for some positive 
constant M 2 ) , 


Propos Itlon 3.3.1. 

If (Hj_) and (H 2 ) are satisfied then all the solutions 
of ( 3 . 1 . 1 ) are bounded. 


The proof Is direct and follows from (3.2.3). 
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Remark 3,3*1 , 

The hypotheses (l^) and (H 2 ) in Lemma 3,3.1 can be 
bartered for the following conditions : 

( ) : r( t, O) is bounded, 

t s+1 - / 

(C*) : Sup / (/ I r(t,T) 1 ds < 

0<t<oo -1 s 

(r(t, s) is defined to be zero outside the region 0<s<t<o>,) 


s+1 


1 . 1 


(C*) : Sup (/ if(T)ii" dt) < 00 where i + i = 1. 


Skco s 


p q 


Indeed, it follows farom (3,2.3) and Holder's inequality 

that 

t 

lx(t)l < fr(t,0)i IXqI + / 1 rCt,s) f(s) Ids 

~ o 

t r 

<M + / / |r(t,T)i if(T)jdsdT, 

“ o t-l 

where M is some constant 

t s+1 

<M + / / lr(t,T:)i if(T:)|dTds 

-Is 

< M+/ ( /"** 1 r(t,T) I / I f (t) i ^ dT:)^/P ds 

“-Is s 

1 

s+1 , / t s+1 — 

< M+ Sup (/ 1 fCx) I^dT)^'^^ Sup / (/ 1 r(t,x) I^)‘5d 

“ 0<s<«» s t>0 -1 s 


Remark 3,3*2* 


It is clear that (H 2 ) implies (C^) but the converse 


is not true. Moreover CC^) and (^ 2 ) hold for r(t,s) = e 


•(t-s) 
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Proposition 3.3.2. 

If (H 2 ) and the hypothesis 

(Hj^) : c(t) IS bounded (i.e. lc(x:)l < for some 
positive constant 

are satisfied, then all the solutions of (3.2.1) are bounded. 
Proposition 3.3.3, 

If (H 2 ) and the hypothesis 

00 

(H^) : Sup / lr(t,s)|dt < 

0<s<<» s 

are satisfied then all the solutions of (3.1.1) are in L 

Proof. From (3.2.3), (H 2 ) and (H^) we have for any T > O, 

T T T t 

; i x(t) 1 dt < / I r(t,0) I ix^ldt + / / |r(t,u)i 1 f (u) i du dt 
o ~ o 00 

00 T T 

< 13^1/ I r(t/0)ldt4-/ (/ |r(t,u) Idt) |f(u)ldu 

o o u 

< 

Remark 3.3.3. 

The hypotheses (H 2 ) and (Hg) in proposition 3.3,3 can 
be replaced by 

(C^) : is [ 0 ,«>) and 

00 

(Cc) : Sup / lr(t,u)! X(t)dt < ~ 
u>0 u 

where X(t) is any nonincreasing function greater than zero. 



36 


Remark 3, 3,4, 

In view of remark 3.3,1 the integrabilxty condition on 

f(t) may be relaxed to Sup / lf(s)l^ ds < «» in Proposition 

u >0 u 

3.3.3, if r(t,s) satisfies (i) r(t, 0 ) is and (ii) 

^ 11 -KL <30 

/ (/ (/ I r( t, s) I dt) du < oo where 1+1 = 1 

-1 u s p q 

instead of (H^) . 

Remark 3.3.5. 

If r(t,0) is integrable over [ 0 ,=°) and if 

oo s 

(Cg) : / / I r( s,T ) f(T) I dx ds is finite 

o o 

then all the solutions of (3,1.1) are integrable over [ 0 /C») , 
For example, if r(t/s) = e ^ ^ and f(t) = e^ then (Cg) holds 
but (H 2 ) does not hold. 

Proposition 3.3,4. 

If the hypothesis (H 2 ) and 

(H 3 ) : c(t) e [o,~) 

are satisfied, then all the solutions of (3.2.1) are integrable 
over [ 0 ,<») • 

Theorem 3.3.1 

If the hypotheses (H^) , (H^) , (H^) and 

(H^) ; The function S>(s) defined by 

ll(s) = ess Sup ia(s,x)l is such that 
0 <T<s’ 

00 — 

/ Jl(s)ds < where is seme constant 
o ~ 
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are satisfied, then the set {0} is eventually unifomly stable 
for (3.1.1). 

Proof. Consider, for O < t^ < t, 

t t 

/ r(t,s)x'(s)ds = x(t)-r(t,t )x(t )-/ [r(t,s)] x(s)ds 

t o o as *- J 

o 

Now (3.1,1) and the change of order of integration yield 


^o ^ 


t t 


/ (/ r( t,s)a(s,T)ds)x(T)dT +/ ( / r( t, s) a( s,t) ds)x( T)dT 


o t 


^o ^ 


t t 

+ / r(t,s)f(s)ds+r(t,tQ)x(t^) =x(t)~/ ^ [r(t, s)] x(s)ds . 


Invoking (3.2.2), we get 


t o t 

x(t) =s r( t,t )x(t ) +/ r(t, s) f(s)ds+J (/ r(t,s)a(s,T:)ds)x(T)dT 

o o 


(3.3.1) 


Prom (H 2 ),(H 2 ) and Proposition 3.3.3 we have 


/ I x(T)jdT < Mg where Mg is some positive constant. 

For a given e > 0, in view of (H 2 ) and (H^) choose 
% ^ ° ^2 ^ '^2^' 2t^Mg ^ ^ ° respectively, so 

large such that 


and 


/ 1 f(s) Ids < for all t > T. 

4Mi 


/ ]i(s)ds < for all t > 7!2. 


(3.3.2) 


(3.3.3) 
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Define a(e) = max { 

Then, for t. > > a(e) > 0 and !x(t^) 1 < 6(e), 

(3*3.2) and (3,3.3) we have 


- ^ 

from (3.3.1) 


, 


t 

!x(t)i < lx(t^)l + / lr(t,s)l lf(s)lds 


+ J° (/ !r(t,s)| ia(s,T)lds) ix(T)!dx 
° ■^o 


< Ml 6 (e) + = e . 

Thus the set {0} is eventually uniformly stable for (3.1*1). 


Remark 3.3.6. 


It can be easily seen that if (Hj_) , (H 2 ) and 

s t 

(C-j) : for s < t, u(t, s)=/ / I r(t, 2’)a(^',T) jdA dx -* O as 

o s 

S *♦ a> 

are satisfied, then the set {0} is eventually uniformly stable. 
In fact assertion still holds if v/e replace (C^) by 

t 00 

v(t) = / / 1 a( s,t) 1 ds dt -♦ O as t -*• «>. 

o t 

3.4 Eventual Uniform Asymptotic Stability 

•In this section the eventual uniform asymptotic 
stability property of the set {0} for both convolution and 
non-convolution equations is discussed. 

Theorem 3.4.1 . 

If the hypotheses (Hj^), (^2^^ ^^3^' ^^4^ 
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(Hg) : r(t,s) •* O along the lines parallel to t-axis 
hold, then the set £0} is eventually uniformly asymptotically 
s-cable for (3«1,1). 


Proof. Since all the conditions of theorem 3.3.1 hold, the 
set £0} IS eventually uniformly stable for (3.1.1). Therefore, 
in view of definition 3.2.3 it is sufficient to show the 
eventual uniform attracting property of £0}. From the 
eventual uniform stability of the set £0}, we have 


for a given T] > 0 there exists 6(^) ^ O and a(T)) > O such thac 
lx(t)l <77 for all t > t^ > a(77) (3.4.1) 

whenever ix(tj^)l < 6(77). 

Prom (H^), for a given > 0 there exists ( |^” ) > O 


such that 


; |f(s)|ds < for all t > Tj^( 

6(77) 

8 Mi 


-‘“I V owr V 


6 (77) 
8M^ 


). (3.4.2) 


0<s<Ti(i^) 


r 6(77) . 

'8nj_ 


Let r(t,s^) rr max _ r(t,s) for all t > 


In view of (H^), for a given 


> O there exists To( i4^ ~) > O 


8M. 


such that 


ir(t,s„)l < IPI for all t > 1,(1^) 


8M^ ” - 3^8M, 


"3'8M, 


(3.4.3) 


Using ( 3 . 4 , 2 ) and (3,4.3), for t > ^ 

get 


r6(77) 
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/ I r(t^s) 1 1 f(s) Ids » / 


T 

'■QM / 




I r(t,s) 1 I f (s) I ds 


+ / 




1 r( t,s) 1 1 fCs) i ds 


< 

- SKj 

Similarly, we obtain 


f/I A ^ ^ ^ ^ M — ^ 

SM^ ^2 + 8M^ • 


(3.4.4) 


/ I r(t,s)i £(s)ds < Ur^ for all t>No=T^(ln^)+To(Tn^) 


2M, 


._..2 — 3 V 4 m ^ M 5 ^ " 4M^M2 ' 


(3.4.5) 

where T 3 (g) is obtained from (H^) for any given g > O. 

Further Proposition 3,3.1 asserts that there exists a constant 
p > O such that lx(t) 1 £ p for all t > O, In view of 
eventual uniform stability of {O}, for P > 0 there exists 
6(P> > O and a(P) > 0. Define 6^ = 5(P) and = a(p). 

Let O < T? < p . Choose T(r)) = N3^+N2-Hi(77) + ^^(l^) and let 

t^ - t^+TCr)), where to ^ cCq > 

Prom (3.3.1) we get 

x(tj^) = r( , t^) x( +/ r(tj, /s)£(s)ds 

t-l 
o 1 


+ S (J r(tj^,s)a(s,T:)ds)x(i:) dx 


(3.4.6) 


o t 


Thus for lx(tQ)l < 6 q# (3.4.6) together with Proposition 3.3.3^ 


(H^), (^^/ (3.4.4) and (3.4.5) gives 

|x(ti)l < 6o + “4— + 25^- ^ 
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Since > a(T)) ^ it follows from (3.4.1) that 


!x(t) 1 <71 for all t > t^ = t^+TC??) , t^ > 


®o* 


Therefore the set {0} is eventually uniformly attracting for 
(3.1.1) and this proves the theorem. 

Theorem 3,4.2. 

If the hypotheses (i^) / (^ 2 )/ (H^) and 

(H^) 1 b(t) = bj^(t)+b 2 (t) where bj^(t) is bounded 

and b 2 (t) is integrable over [O/") 

are satisfied^ then the set £0} is eventually unifoinnly 
asymptotically stable for (3,2.1) • 

Proof. As in the proof of theorem 3.3.1, from (3.2.1), (3.2.4. 
and for all t > t^, we have 

t t t 

y(t) = c( t-s) f( s) ds+/°( / c( t-s) b( s-T)ds)y('c) dx 

t o t 

o o 


(3.4.7) 


Consider the integral 

t X. 

o t 

1/ (/ c( t-s) b( s-T)ds)y(x) dti 


o t. 


4- 

t o 

< / (/ lb(s-T)l iy(T)ldT) |c(t“s)lds 

” t o 
o 

t s 

< / (/ |b(s-T)| [y(T)id't) ic(t-s)lds 
“ -^o ° 

t t— s 

< / (/ tb^(t-s-T)l ly(T)idT) ic(s)l ds 


to o 


+ / ^(s) 1 c(t-s) J ds 

o 


(3.4.8) 
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s 

where l(s) =/ lb,(s-T)| ly(T)ldT. 

o 

From (H^) / Proposition 3.3,4 and (H^) / we have for a given 
e > 0 there exists a T^Ce) > o such that 


t t-s 

/ ( / I bj_( t-S“T) 1 


i y(T ) 1 dx) I c( s) i ds < ^ for t ^ ^ ^3 

(3.4,9) 


In view of hypotheses and equation (3,2.4) it is clear that 

c'(t) is "bounded and hence c(t) is uniformly continuous. 

Also by (H^) and Proposition 3.3.4, f(s) e L [o,«>) and is 

positive. Since c(t) is uniformly continuous and 
1 ~ 

f(t) is L [0,«>), / W(s) lc(t-s)ids is uniformly continuous 

o ^ 

a'nd moreover it is integrable because lS(s) is L [o,“>) and 

(H^). Now the application of a lemma of Barbalat L*^f. 13, 

OO 

p. 89j yields that S ^(s) lc(t-s)| ds -• 0 as t -* «»• 

o 

Therefore for a given e > O there exists a T 2 (e) > O such that 


I 

o 

By (H 2 ) 
that 


iHs) lc(t-s) Ids < !■ for t > t^ , 
for a given e > 0, there exists 



T = (-^) 

4M^ 


(3.4.10) 
> O such 


/ 1 f ( s) 1 ds < for all t > T, • 

4i\ 

Define a(e) = max {T 2 ,T 2 ,Tj^(-;; 7 -) } and = 

4M^ 

Then for t > "to - ^ ^ iy(tQ)f < 6(e) 

( 3 . 4 . 7 ) - (3.4.11) give 


(3.4.11) 



jr ecjiations 
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{y(t) I < 6(e) + + I + I = e. 

4M^ 

Hence the set {0} is eventually uniformly stable for (3.2.1). 

Since c(t) is uniformly continuous and is integrable, c(t) 
as t -*■ <», 

Now proceeding as in the proof of Theorem 3.4.1, it 
can be easily seen that the set {0} is eventually uniformly 
asymptotically stable for (3.2.1). 

3.5 Perturbed Equations 

The eventual uniform asymptotic stability of the 
set {0} for (3.1.2) is discussed in this section under the 
class of peirterbations that satisfy the following hypothesis 

(Hg) : lg(t, z)l < \(t) (1+lzl) for t > 0, 1 zl < » 

where X(t) s [0/«») . 

Theorem 3.5.1. 

Let the hypothesis of Theoron 3-4.1 together with 
(Hg) be satisfied. Then the set {0} is eventually uniformly 
asymptotically stable for (3.1.2). 

Proof. Hypotheses (H^) and (H^) imply that 

ifj rCt,s)l < Mg C3.5.1) 

where Mg is scatie positive constants 

Prom the variation of constants formula, the solution z(t) 
of (3.1.2) can be expressed as 
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2(t) = x(t)+r(t,tQ) ( z(tQ)--x(t:Q) ) 

■t 

o t 

f I (J r(-C/U) aCu^T ) c3u) g(T#zCT:))dT 
■t ^ 

- / --r r(t|.T:) z(T))dT 

4- 3^ 

t t 

+ / °( / r( t,u) aCu/";) du) ( z('c)-x(t ) ) dT, (3,5.2) 

o t 

o 

where is some arbitrary point and x(t) is a solution of 
the linear equation (3,1.1), For = O ^nd for z(o) = x(0), 
the equation (3.5.2) takes the form 


t 

2(t) = x(t) - / 
o 

then the hypothesis 


L-. 

dT 


(Hg) 


r(t,T) g(T, z(T)dT , 

together with (3.5,1) 


gives 


t 

lz(t)l < M^4Mg / X(t) 1z(t)i dT 
where is a positive constant. 

Since 5v.(t) e L'*' [O^o*) , the application of Gronwal 1-Bellman 
inequality yields that every solution z(t) of (3,1,2) is 
bounded. If x(t) is a solution of (3.1.1) such that 
x(tQ) = z(tQ) then (3.5.2) together with (3.5.1)# (H^^)# (H^) 
and (Hg) gives 


t 

lz(t)-x(t)l < p(t#tQ) + q(t#t^) / Iz(t:)-x(t)i 

t t 

where p(t#t ) = Mq(/ i.(u)du) + Mg / X(u) du 

t 

o o 


dx 

(3.5.3) 


and 
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t 

q(tyt^^) =M 3 ^(/ £(u)du), Mg^Mg being positive constant":, 

"^o 

T+- be c 2 early seen that pCo/t^) and q(t, t^) are bounded 

and increasing in t« Further/ pCt,t^) /q( t/t^) -» O as t^ -* "• 

Applying Gronwall- Bellman type inequality [cf. 53/p.l6] to 

(3.5.3) we obtain ^ 

^ S q(s/t^)ds 

lz(t)-x(t)| < p(t/t )+q(t,t ) (/ p(S/tQ)ds)e° 

< p('C/tQ)+q(T,tQ)Nj^T e ^ for t e [o "c), 

where Nj^/N 2 are bounds on pCt/t^) and q(t/t^) respectively. 

NT 

Let = pCT/t^) +q(T /t^) e ^ . Then, it is clear 

that - 0 as t^ 

Thus the hypotheses of Lemma 3.2.2 hold with K(.) = O and 
hence the set {0} is EvUAS for (3.1.2) • 



CHAPTER - 4 


LIAPUNOV- RAZUMIKHIN TECHNIQUE 
4 •! Introduction 

Liapunov second method or the direct method involving 
an energy like function has come to stay as a powerful tool 
in the qualitative study of ordinary differential equations. 
Over the years^ this method has also been extended to 
functional differential equations (FDE) by various authors. 
While Liapunov functions are employed in the study of ordinary 
differential equations/ more generally Liapunov functionals are 
used in studying PDE (cf, [l 9]/[5 4 ])* As Vol terra integro— 
differential equations (VIDE) can also be treated as PDE, 
such a study can be carried over to these equations. Further 
Liapunov fxmctionals have been constructed exclusively for VIDE 
by Burton ([l] / [ 2 ],) in order to study the stability and 
uniform stability properties of VIDE, As distinguished frran 
this line, Rafumikhin [ 41 ] obtained stability properties of 
FDE using Liapunov functions rather than functionals, Ihis 
technique of Rammikhin involves an estimation of the derivative 
of the Liapunov function along the solutions of the system of 
equations on a certain set in which the solutions satisfy a 
specified inequality. In recent y^rs a good number of papers 
Ccf* £153^ [21]# [243# [473, [493 ) have appeared employing this 
technique and various results have been accumulated. In 


47 


particular, Seifert [47], [48] and Grimmer and Seifert [15] 
have studied the properties of solutions of VIDE such as 
stability, boiandedness, asymptotic stability, uniform 
stability and so on, by studying the corresponding properties 
of FDE. In this chapter Liapunov- RaxumiMiin technigue is 
further explored to yield L —stability of the zero solution 
of VIDE, This has been achieved as an application of a more 
general result on PDE which has been derived in section 3, 

Section 2 deals with the preliminary results that 
are necessary for proving the mam results. For completeness 
some of the results of Seifert are included in this section, 

4.2 Preliminaries and basic results 

In this section the etiployed notations are eaqilained 
and then it is indicated how a VIDE can be thought of as an 
ordinary differential equation involving an interval of delay 
which becomes unbounded as t -* +<»« In addition two theorems 
of Seifert [47] are given \diich ensure the stability of the 
point X = O and boundedness of the solutions,for (4,2,1) , 

For an element x of n”, by | x| we mean the usual 
euclidian norm# We denote ty a function continuous on 

the interval 0 < s < t to IStP’ and by S^ we mean the set 
£Xj.(,)} of all such functions. If x{s) is a function defined 

XX 

and continuous onO<s<»to H, then for each fixed -t,. 

O < t < «>, this function defines a member of S^ given 

by x(s), 0 < s < t. We call this function x^(,) a segment of 
x(s) • 
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For fixed t > O, let F(t,Xj-(,)) be a function on 
to • For each function x( s) continuous on 0 < s < «> to 
we assiime that F(t^x^(,)) is continuous in t/ where 
x^( •) is a segment of x(s). 

By a solution of the equation 

x'(t) = F(t,x^.(.)) (4,2.1) 

we mean a continuously differentiable function x(s) on 
O < s < <» such that (4,2.1) is satisfied on 0 < t < «> for 
Xt.(.) a segment of x(s). Let x{t^XQ) be any solution of 
(4.2,1) existing for all t > O such that x(0,x^) = x^. 

Definition 4,2,1. 

The point x = O for (4.2.1) is stable if given 6 > O, 
there exists a 6 (s) > O such that if 1x^1 < 6 (e)/ then every 
solution x(t/XQ) of ( 4 , 2 , 1 ) is defined for t ^ 0 and satisfies 
lx(t/>:^) 1 < e for t: > O. 

Definition 4.2.2, 

The point x = O is L^-stable, 0 < p < », for (4,2 •!) 
if it is stable and if there exists a 5 ^ > 0 such that 

00 

/ |x(t/ 3 ^)l^ dt < <» whenever Ix^^l < 6 q* 

o 

We shall now consider the VIDE of the form 

t 

x'(t) = H(x(t)) + / G(t,s^x(s) )dS/ x(0) = 3 ^ 

(4.2.2) 

for t > O 0 Here H(x) is continuous on 3R , and G(t/S/x) is 
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continuous on x 

It can benoted that if x(t) is a solution of (4.2.2), then 
it IS also a solution of (4.2.1) v^diere 

t 

I’(t,g(,)) =H(§(t)) + / G( t/S/gt's) )ds 

o 

and x^(.) IS a segment of x(s), 0 < s < <». 

We shall now state the following results due to 
Seifert [47] which are useful in our subsequent discussion. 

Theorem 4.2.1. 

Let there exists functions u(s), v(s), and f(s) 
continuous for s > O and such that u(0) = v(0) = 0 , u(s) 
is increasing, f(s) > s for s > 0, and suppose V(t,x) is a 
real- valued function continuous in (t,x) for t > 0 and x in 
D, an open subset of 3R containing the zero vector. Let V 
satisfy : 

( i) u( I xl ) < V( t/x) < v( 1 xi ) for t > 0, X in Dj 

( ii) v'(t/x(t)) < O for any solution x(t) of (4.2.1) for 

which x(s) is in D and f (V( t,x( t) ) ) > V(s,x(s)) for O £ s < t. 
Then the point x = O is stable for (4.2.1). 

Remark 4.2.1. 

In Theorem 4.2.1 and onwards, by V'(t,x(t)) we mean 
[V( t+h , x( t+h) ) - V(t,x(t))] 


lim Sup 
h-O**- 


h 
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Theorem 4, 2 •2* 

Let there exist a function V satisfying the hypotheses 
of Theorem 4*2 *1 except that now D = IR^ and u(r) -► oo as r ■* «>. 
Then the solutions of (4.2.1) are bounded. 

4.3 L^-Stability of Functional Differential Systems 

In this section# we shall first obtain a general result 

yielding L^-stability of the point zero for the FDE (4.2.1). 

2 

Then this result has been applied to derive L -stability 
properties of the zero solution of an integro-differential 
system in Section 4.4. 

Theor^n 4,3.1. 

Suppose V : HR*" xIr’^ 31^ is a continuous function and 
satisfy the following conditions : 

(i) u(lxl) < V(t,x) < v(lxl) for all (t,x) 

where U/V : HR* ■* IR"^ are continuous functions such that 
u(0) = O = v(0) and u(s) is increasing to as s increases 
to 

( ii) For any solution x(t) of (4,2.1) 

V'(t#x(t)) < - C|x(t)|^ # p > O# C > O# 

whenever f(V(t#x(t))) > V(s#x(s)) for s in any neighbourhood 
(t>-n#t)# 1 ) > 0 , and where f(s) is any continuous function for 
s > 0 such that f(s) > s for s > O, 

Then the point x = 0 is L^-stable for (4.2.1). 
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Proof. Define 

t 

r(t) = V(t,x(t)) + C / lx(t,x )l^ dt 

o ° 

Clearly r(0) = V(0 ,Xq) and r(t) is positive, 

We claim that r(t) < r(0) for t > 0, Suppose not. Then 
there exists a tj^ > O such that r(^) > r(0). Since r(t) is 
a continuous function of t^ there exists a > O and a t 2 > t 
such that r(t) < r(0) on , r(tj^) = r{0) and 

r(t) > r(0) on (t^,t 2 ). Therefore for some t* e 
there exists an T) > o such that 

r'(s) > 0 for s e (t*-h,t] . (4,3.1) 

Further for t*- 77< s < # # 

r(s) < r(t*). (4^3,2) 

From hypotheses and Theorem 4.2.2 it follows that 

Sup IxCt/X^)! < 

CK 

where Mj^(xj^) is a positive constant and x(0,Xq) = x^. Frcxn 
continuous dependence of solutions on the initial values 
we get < M for 13^1 < 1/ where M is a constant 

independent of Xq. 

Hence the monotonocity of v yields 

v(|x(t^3^)|) < v(M)^ for t > O, JXqI <1. 

. 83lt7': 
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Let e be given* We may assume that e so small that u(fi) v(m) « 
Then there is a positive number a = a(e) such that 


f(s) > s+a for s e [u(e),v(M)] • 


(4.3.3) 


From the hypotheses and Theorem 4.2.1 it is clear that 
the point x = 0 as stable for (4,2.1). (4.3,4) 

Therefore^ for a gaven L = there exists a ^ > 0 

such that 


Ix^^l < implies |x(t^3^)| < t > O, 

In particiilar for s e (t^-Tj^t^') and for 1x^1 ^ 6 q = 
we have 

Jx(s,x^)lP < 2!^ • 

This implies that 

C / Ix(t/Xq) dr < ^0 (t*-s) 

3 

< I < a (4.3.5) 

From ( 4 . 3 . 2 ) and the definition of r(t) , we get 

V( t*/x( f^") ) + C J 1 x(t/X )l^ dx > V(sjrx(s)) ( 4 . 3 . 6 ) 

s ° 

For a given 6 > 0 there exists a ~ 6j_(®) ^ O such that 
v(e) < u(6i). Since x = O of 4,3.1 is stable, it is clear 
that u(6i) 1 V(t,x(t)) < v(M), Hence from (4.3.3), we have 
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f(V(t*^x(t--) )) > v(t*^xCt*)) + a 

> v(t*,x(t*)) -rC / |x(T,x)l^dT: 

s ° 

( f rcsn 4,3.5) 

> V{s,x(s)) (from 4.3.6) 

rpT.^ 

s condition (ii) of hypotheses gives 
V'(tOx(t*‘-)) < - C!x{t*)|P 

which in turn implies r'(t*) < O contradicting (4.3.1). 
^^refore, there exists no > O such that r(^) > r(0), 
^®hGe rCt) < r(0) for all t > O. Thus we have 

t 

^(t,x(t)) + C / 1 x(t#x^)|^ dT < V(0,x^) and therefore 

o 

t n 1 

/ ix(T/X^)r dT < ^ V(0 ,Xq) for ix^^l < 6^. 

'^his together with (4.3.4) proves the theorem, 

Integro~Differential Systems 

As an application of Theorem 4.3.1, we consider an 
^^^®gro-differential system 

t 

x'(t) = Ax( t)+h(t,x(t) )+/ g(t,s,x(s) )ds, x{0) = 

o 

(4.4.1) 

A is a real nxn constant matrix, h is continuous on 
^ satisfying 


lh(t,x)l < w IxU 


(4.4.2) 
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and g(t,s^x) is continuous on 3 R'*' xIr"^ x for t > s > 0^ 
and satisfies 


lg(t,s,x)l < L(t,x) ixl (4.4.3) 

t 

with / L{t,s)ds -* 0 as t - oo. 
o 

Theorem 4.4.1. 


Asstime that 

( i) A IS a stable matrix^ 

(ii) (4,4.2) and (4.4,3) hold. 

2 

Then for sufficiently small , x = 0 of (4.4.1) is L —stable. 

Proof. Since A is stable, it is clear (cf. [44]) that there 
exists a positive definite symmetric matrix B such that 

a'^B + BA = -I 


T 

where A is the transpose of A and I is the identity matrix. 
T 

Define V = x Bx. Let a and 0 be respectively the smallest 
and largest eigenvalues of B. Then we have 


'> 0 T 9 7. 

or 1 xr < X B X < r i xr 


(4.4.4) 


for X e IR^, Hence the condition (i) of Theorem 4.3.1 is 
verified. To show that the condition (ii) of Theorem 4.3.1 
holds for jj. sufficiently small, let 


a 


^ ^ 2lBi$ 


n 

S Ibiil* ® 

i<j=l ^ 


(4.4.5) 


(by) 


where i Bt 
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Choose f(s) = q^s where q > 1 is such that 


^ 2ql Bl 0 (4.4.6) 

Then for any rj > o and any solution x(t) of (4,4.1) such that 


X (s) Bx(s) < q x^(t) Bx(t) for s e [t-r),t] , (4.4.7) 

we have 


2 2 2 2 2 

q 3 !x(t)| > a |x(s)i"^ for all t > T) . (4.4.8) 


For a solution x(t) of (4,4.1) with V(t,x(t)) = x^(t)Bx(t), 
it follows for t > 0 that 


V( t,x(t) ) < -.ix(t) + 2|B| ix(t)l |h(t,x(t))i 

t 

+ 2|Bl |x(t)| / t g(t,S/x( s) ) I ds, (4.4.9) 

o 

Fix jLij^ > 0 such that 


(2IB^, UJ. < 1. (4.4.10) 


In view of (4,4,3) there exists a rj > O such that for t > TJ , 
t 

/ L ( t ^ x) ds ^ m 
o 

Since 

t t 

/ I g(t/S^x(s) ) 1 ds < / L(tjrx) Ix(s)|ds 
t-v “ t-n 

t 

< Sup x(e) / L(t^s)ds 
. “"eert-n^ti o 

■* (4.4.12) 
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Prom ( 4 . 4 , 11 ) and (4.4.12)/ for a solution x(t) of (4.4*1) 
satisfying (4.4,8), we have 

^ 2qBjl^ I Bl 

21B1 lx(t)l S lg(t,S/x(s))lds < i |x(t)i^ (4.4.13) 

t-f) ■” ^ 

Further (4.4.3) gives 

t— T) ^ 

j. i g(t/S,x(s) )i ds < M / L(t/s)ds - O as t « (4.4.14) 

o o 

Thus ( 4 . 4 . 9 )/ ( 4 . 4 , 10 )/ (4.4.13) and (4.4.14) lead to assumption 
( ii) of Theorem 4.3.1. Hence the application of Theoren 4*3.1 
yields the desired result. 


CHAPTER - 5 


ANALYSIS OF LINEAR SYSTEMS 

5 *1 Introduction 

The study of linear systCTis is vital in order to 
penetrate into the complexity of non-linear systems. Besides^ 
the linear systems themselves are very often fitted as 
convenient mathematical models for various physical phenomena. 
Thus it is Important to accumulate more information regarding 
the behavioural pattern of linear systems. This chapter has 
been exclusively devoted for the analysis of linear systems. 

It is well known that, a linear autonomous ordinary 
differential system is asymptotically stable if all the 
eigenvalues of the corresponding coefficient matrix have 
negative real parts Ccf. [44, Chapter 3]). For non-autonomous 
systems, with an addition of Lipschitz condition on the 
coefficient matrix, similar results have been e^qpounded in 
[ 4 ], [ 5 ], Thus v^ile studying Vol terra integitJ-differential 
systems (VIDS},be it through Lyapunov second method (cf. [l]/ 
[2], [471/ [48]} perturbation analysis (cf. [42], [44]) 

it has Invariably been assumed that the coefficient matrix 
is stable. But this amounts to requiring the stability of 
the corresponding ordinary differential equation. Notable 
exceptions that have dispensed with the stability condition 
on the coefficient matrix have been the works of Levin [25] / 


* 
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Grossman and Miller [l7]j> Griiraner and Seifert [iS] , and 
Burton [sj. In [25] this has been done by defining a suitable 
energy function vdiile in [17] the integrability of the 
resolvent function of VIDS has been characterised by a 
transformation condition similar to that given in [39] for 
Volterra integral equations. In [l5] / the same has been 
achieved by studying the resolvent of a transformed equation. 
Quite recently in [ 3 ] the conditions involving the anti- 
derivatives of the kernel are assumed. Motivated by the 
interesting nature of this problem and the above cited works,^ 
an attempt has been made in Section 2 to study the asymptotic 
behaviour of solutions of VIDS -when the coefficient matrix 
is not necessarily stable. Our approach has been through 
an equivalent VIDS Involving an arbitrary function satisfying 
certain assumptions and thus paving the way for the new 
coefficient matrix (corresponding to equivalent VIDS) to be 
stable. 

The results obtained in Section 2 have been used to 
study a linear ordinary differential system and the conditions 
for integrability of its solutions vhen they are bounded and 
vice versa are discussed in Section 3. 

A line of thought in recent years has been that the 
dichotcanies rather than Liapunov characteristic expon^ts form 
a proper base for the study of asymptotic behaviour of solutions 
of non-autonomous differential systems. Massera and 
Schaffer [32], Coppel [5*6,7] * Saker and Sell [46] among 
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others have discussed dichotomy propeirties of ordinary 
differential systems and several results have been reported. 
These properties have been extended for Vol terra integro- 
differential systems in Section 4. 

5*2 Stability properties of solutions 

In this section it is proposed to study the stability 
properties^ more specifically the exponential stability/ of 
the solution x(t) of 

t 

x'(t) as A(t)x(t) + J K(t/S)x(s)ds + F(t)/ (5.2,1) 

o 

wherein A(t) and K(t/s) are nxn matrices defined and continuous 
over O < t < w and o < s < t < *» respectively and xCt) and 
F(t) are n-vectors with F(t) continuous over 0 £ t ^ when 
the matrix A is not necessarily stable. Our main approach 
is by way of an equivalence theorem (Lemma 5 •2*1) that has 
the potential to supply us with a stable matrix B correspondin’^ 
to A, The efficiency of this approach has been dsoonstrated 
through a suitable ejciraple. 

The following lemmas are useful in our subsequent 
discussion, 

Lonma 5,2*1 

If $(t/S) is any continuously differentiable function over 
O < s < t <»/ then equation (5.2,1) is equivalent to 

t 

x'(t) » B(t)x(t) + / L(t/S) x(s)ds + H(t) 

o 


(5.2.2) 
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where B(t) = A{t) - 

L(t,s) = K(t,s) + $g(t,s) + f(t,s)A(s) 

t (5.2.3) 

+ / §(t,u)K(u,s)du, 

s 

and 

H(t) := F(t) + #(t^0)xQ + §(t,s)p(s)ds. 

o 

Proof. Let x(t) be any solution of (5.2.1) existing on the 
interval 0 < t < <». Consider the identity 

^ t 

/ fs(t,s)x(s)ds = #(t,t)x(t)-§(t,0)x^ - / $(t^s)x''(s)ds. 

° o 

Substituting for x^(t) from (5.2.1) and using Pubinis theor^^ 
we get 

^ t 

/ #s^^-'S5x^s)ds = §(t,t)x(t) - §(t^0)x -/ $(t,s)F(s)ds 

° ° o 

. (5.2,4) 

■t t t 

- / §(t,s)A(s)x(s)ds - / (/ ${t,s)K(s,T)ds)x(T)d'r 
o o T 

Therefore, it follows fran (5.2.1) - (5.2.4) that 

t t t 

/ L(t,s)x(s)ds = / K(t,s)x(s)ds + / if_(t,s)x(s)ds 


■t t t 

+ / f(t,s)A(s)x(s)ds + / (/ ${t,u)K(u,s)du)x(s)ds 
o OS 

= x'(t) - A(t)x(t)-F(t)+f(t,t)x(t)-§(t,0)3^ 
t 

- / $(t,s)F(s)ds 
o 

= x'(t) - BCt)x(t) + H(t) and x(o) = 
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Thus every solution of (5.2.1) is also a solution of (5.2.2). 
Conversely, if y(t) is any solution of (5.2.2) then by taking 
$(t,s) = 0 we have from (5.2.2) and (5.2.3), L(t,s) = K(t,s) 
and H(t) = P(t) and hence y(t) satisfies (5.2.1). 

Remark 5.2.1. 

It can be readyly observed that when $(t, s) is the 
differential resolvent corresponding to the kernel K(t, s), 
then the equation (5*2.2) together with (5.2.3) gives the 
usual variation of constants formula (See Grossman and 
Miller [16]). 

Lanma 5.2.2. 

If B(t) is any nxn matrix function which commutes with 

its integral and satisfies the condition 
t 

/ B(T)dT 

I gS I 1 M O < s < t, (5.2.5) 

then the solution x(t) of (5.2.2) with x(0) =3!^ satisfies 

|x(t)| <M{x^| + M / 1 H(T)| dT 

o 

tM / (/ lL(T,s)j d-C) lx(s)! ds. (5,2.6) 

o s 

Proof. Multiplying both sides of equation (5.2.2) by 
t 

exp (— / B(T)dT) and rearranging the terms, we obtain 
o 

t t 

-/ B(T)d^ -/ B(T)dT t 

(e ® x(t))* = e ° [H(t)+/ L(t,s)x(s)ds3 • 

o 
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(v) l$g(t^s)| < N e Ny 6 are positive constants 

(vi) for this $, the hypotheses of Lenma 5.2.2 hold, 

and 

( vii) y ^ p > a, 6><i and < ct where 



Then every solution x(t) of (5.2.1) tends to zero exponentially 

as t -* CO, 

Proof. 

In view of lemma 5.2.1 and for a function § satisfying 
the conditions (iv), (v) and (vi), it is enough if we show 
that every solution of the e<^ation (5.2.2) tends to zero 
exponentially as t -► «». Since P(t) = O, equations (5.2.3), 
(5.2.6) and condition (iv) imply that 


lx(t)| < MI 3^1 +ml13^I / dr 

o 


t t 


at 


+ M / (/ e |K(t,s)+« (x,s) + $(t,s)a(s) 
os ® 


+ / f ( t,u) K(u,s)dujdt) ix(s)lds 
s 

<M,v - Ti - 


t t 


OCT. 


+ M / (/ e iK(r^s)i dr) lx(s)i ds 
o s 

t t 

+ M / (/ e |# 3 (T^s)ldT) ix(s)|ds 
o s 



64 


t t 


+ M J (J ^ !f(T,s)i !A(s)l dT) tx(s)|ds 
o s 

t t "C 

+ M/ [ / e®^{/ I §(t,u)i 1K(u,s) idu)dTl |xCs)}ds. 

OSS 


Further, from (ia) and (vii) we have 


e“^ix(t)! <M15^1 + + M /c/V^‘^-’®^K(T:,s5jdT)e“^|x(s) d£ 

O 55 


t t 


o s 


t t 


+ M / (/ IS,(T,s)ldT) e®® lx(s)!ds 


o s 


t t 


+ MA^ / (/ e®^'^“®^$(T,s)idT)e“®|xCs)ids 


o s 


*fc *t 

+ M / [ / I 4('E/u) 1 1 K{u,s) I du)dT j 

s 


o s 

Jx{s)jds« (5-2.7) 

i 

Now we shall estimate each integral on the right side of (5.2.'/)e 
Us \vig condition ( iii) , we get 

/ gttCT-s) 


< j / e-O-aH 

s 


(5.2.8) 


( 8 -a) 

Condition ( ii) gives 

/ \§A’t,s)l dx 

s ® 

< H /e-(6-a)T ^ 

iS 

Cs #2 
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Similarlyy estimating the other integrals by usxng ( iii) 
and ( iv) , we obtain 

(5.2.10) 

and / | $(i: ^u) j i K(u^s) i du) dr 

s s 


< LJ 


LJ 

(y-3)0-a) 


g~(0-a) (t-s) 


. LJ 

(y-0nr-a) 


g-(r-a)(t-s) 


(5.2.11) 


The inequalities (5.2.8) - (5.2.11) and (5.2.7) lead to 


at MLix^l t; mijA 


MLA 


MLJ . MLJ _-(y~a)(t-s)^ as, , 

( 0-a)(y-a) (y-8) (y-a) J e i x( s) i ds 


i»iV h +Tr^] 


MN 


MLA 


MLJ T as . , , 

■*■ r g- S ;Tr- 373 ® IxCsJids 


<M|3!^I [1+Xt^]+^ % e^^ixCs)! 


ds. 


Applying Gronwall-Bellman inequality we obtain 

. at 

e |x(t)l < MIXqI (1 + ' {ylg} ) ® ° • 
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This implies that 

, . T (cc “a)t 

lx(t)| < mi3<^i (1 e ° 

Thus in view of condition (vii), the result follows. 

Remark 5*2.3« 

It is possible to select a matrix function §(t/s) 

satisfying conditions (iv)/ (v) and (vi) of theorem 5.2.I. 

For example if $(t,s) = L I, then N = ^ and 6 = r, 

y 

i(t^t) being a constant matrix in this case/ the estimate 
(5.2*5) is guaranteed if A(t) is a constant matrix and B is 
negative definite. 

Remark 5.2.4 

Basically it is the condition < oc" which controls 

the asymptotic nature of the solution x( t) * A look at the 
composition of reveals that \diile so choosing y and § 
much av;ay from c: we can nullify the effect of the last three 
terms in the first term [ essential term 

with which we have to reckon with* Therefore if ^ is so large 
as to exceed (a +MJ)/a/ then only would be less than a. 

Thus we see that the attenuation required on the kernel 
K(t/s) is linked with the constant a in (5-2.5). This 
conclusion implicitly assumes the estimate (5.2.5) . Such an 
estimate is always possible when the transformed matrix B is 
constant and negative definite. 
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Renark 5«2*5, 

If F(t) is not zero in theoren 5.2.1/ still the 
solutions of (5.2.1) tend to zero provided FCt) is integrable. 
This is an immediate consequence of variation of constants 
formula (cf. [l6]) and theorem 5.2.1. 

Remark 5.2.6 . 

In [ 3 ] , a condition of the type (5.2.5) has been 
def _ 

used for the matrix Q = [A(t) - G(t/t)] / where G(t/s) is 

the antiderivative of the kernel K(t/S)(i,e. = K(t/s}) 

As such Idle matrix B in our study allows more flexibility due 
to the arbitrary character of the function #(t, s). Further/ 
our approach is entirely different and the analysis in [ 3 ] 
can be applied to the equation (5.2.2) in order to obtain 
sharper estimates. Thus our theorem 5.2.1 is in addition to 
the theorem 2 of [3]| rather than a substitute for it. 

Example 5.2.1 

■•loC t*^s) 

In (5.2.1) (scalar case), let K(t/S) = e , 

A(t) = lOe ^ - 5 and F(t) = O. Choose $(t/s) = lOe , 

then a = 5 and 3.25. Thus all the conditions of 

Theorem 5.2.1 are satisfied. 

5.3 An ordinary differential equation 

In this section, an ordinazy differential equation 
is considered as a special case of Vol terra integro-dif f erentie 1 
equation (5.2.1) where the kernel K{t/s) s O. Lenina 5.2.1 has 
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been made use of to study L^-like properties of solutions of 
this equation •when -they are given bounded and vice versa# 

Consider ■tJie equation 

x^(t) = A{t)x(t), x(0) = (5-3.1) 

Invoking Lemma 5.2.1 with K(t/s) = 0 and P(t) = 0^ the system 
(5.3.1) is equivalent to 

t 

X (t) = [ A( ( t, t)l x(t) -f / [ §g( s)A(s) 1 x(s}ds 

o 

+ §(t,0)x^ , 

This in turn reduces to 

t 

#(t,t)x(t) = #(t,0)xQ + / [$g(t^s)+$Ct,s)ACs)] x(s)ds- 

o 

( 5 . 3 . 2 ) 

Further/ if we select the function §(t/s) such that §(t/t) is 
non— singular for each t > 0/ then (5.3.2) -takes the form 

t 

x(t) = G(t) + / U(t/s)x(s)ds/ (5.3.3) 

o 

where G(t) = $'”^(t/t) f(t/0)x^ 

and 

U(t/S) = #'“^(t/t) [$g(t/S) + $(t/s)A(s) 3 • 

Clearly (5.3.3) is in the form of Vol terra integral equation. 
Theorem 5.3.1. 

If the matrix function §{t/s) is such that 
(i) U(t/t) is non-singiilar for ^ch t > 0 
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t 


(ii) Sup / |U(t,s) - U(t,t)l ds < « 

0< tCoo o 


and (iii) Sup lG(t)t < «>, 

then every bounded solution x(t) of (5* 3*1) satisfies 

t 

Sup U(t,t) 1/ xCs)ds| < <» (5«3.4) 

0 <t<<» o 


Proof. From (5.3.3), we have 


t t 

x(t) - U(t,t) / xCs)ds = G(t) + / [u(t,s)-u(t,t)] x(s)ds 


and hence 

t b 

1 x( t) —UC t, t) / x(s)dS| < lG(t)| + / I U( t, S) — U( t, t) 1 lx(s)l C3 
o “ o 

This together with the hypotheses gives the result. 


Theoraa 5.3.2. 

If x(t) is any L^- solution of (5.3.1) and if f(t, s) is 
a matrix function such that 


(i) 

u(t,t) 

is bounded 

(ii) 

Sup 

|U(t,s)-U(t,t) 1 < 


0< s< t<«> 

( iii) 

Sup 
o< t^®° 

jG(t)l < “/ then 


x( t) is bounded. 

The proof is similar to that of theorem 5.3.1 and hence omitte< 
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3 • 3 • X 

For a choice of f(t,s) = e 1/ the conditions of 
Theorem 5.3*1 reduces to 

{ i) (A(t)+l) is nonsingular, 

(ii) Sup 1/ A(s)dsl < “ and 

0<t<00 Q 

(iii) Sup t [A(t)+I] < oo respectively. 

0<t<<» 

For the same choice of $(t,s), the assumptions (i) 
and (ii) of theorem 5.3.2 respectively takes the form 

( i) ( A( t) +I) is bounded and 

(ii) Sup A(s) < ®. 

CKs< t^“ 

Example 5.3.1 

Consider the scalar differential equation 



Take $(t,s) = then all the conditions of 

Theorems 5.3.1 and 5.3.2 are satisfied. The solution x(t) 
of (5.3.5) with x(0) = Xq IS given by x(t) = XQ(l+t)e 
which is integrable and bounded* 
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5.4 Dichotomies of solutions 


While studying the stability and asymptotic stability 

I? 

properties of autonomous linear differential equations, we find 
that these properties are basically characterised by the sign 
of the real parts of the eigen values of the coefficient matrix. 
However an automatic generalisation of this approach does not 
carry over to non— autonomous equations for studying the 
corresponding properties of uniform stability and uniform 
asymptotic stability. An example indicating this failure has 
been given by Hoppensteadt [20], Therefore for such equations 
the study of dichotomies asstmes inportance. 


In this section after giving necessary definitions and 
a lemma, the dichotomy properties of systems of Vol terra integri- 
differential equations are studied. 


The solution x(t) (cf, Grossman and Miller [19]) of the 


system 


x^(t) =s A(t)x(t) + / K{t,s)xCs)ds, x(0) = x^ 

o 

(5.4.1) 


is given by x(t) = RCt,0)x^ 


(5.4.2) 


where the nxn matrix function R{t,s), o < s < t < <», is the 
resolvent corresponding to the functions A(t) and K(t, s). 


Lemma 5.4.1 

R(t,0) is invertible for each t > 0. 
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Proof Suppose not, then a a tj^ > 0 such that R(tj^,0) 
is singular. Therefore the syst^ of equations R(tj^,0)x = O 
has a non-zero vector, say as a solution. If x(0) = x^ is 
any initial point of the solution then by (5.4.2) we have 

xCtj^) = R(t3^,0) Xq (5.4.3) 

Also if X = x^ + X then x / 3 ^ and 

R(t]^,0)x = R( tj^,0)xQ+R( t^,0)x = x(tj^)+0 = x(t 2 ^) (5.4.4) 

Thus (5.4.3) and (5.4.4) indicate violation of uniqueness of 
solutions of (5.2.1) and hence a Contradiction. Thus the 
lemma is proved. 

Definition 5.4.1 

System (5.4.1) is said to possess an exponential 
dichotomy if there exists a projection matrix P and positive 
constants M,N,a, 3 such that the following conditions are 
satisfied. 

jR(t,0) PR"^(s,0)i < M for t > s 

(5.4.5) 

and lR(t,0) (I-P)r’’^(s, 0) I < N for s > t 

Definition 5.4.2 

System (5.4.1) is said to possess bounded growth on 
an interval J =3R if, for seme fixed r > O, there exists a 
constant C > 1 such that every solution x(t) of (5.4.1) 


satisfies 


Lemma 5.4,2. 


System (5.4,1) has bounded growth if and only if there 
exists constants L, r such that its resolvent matrix satisfies 

lR(t,0)R“^(s,0) 1 < L for t > s (5.4.6) 

Proof : Suppose (5.4,1) has bounded growth. Then for 
s < t < s+r 

!R(t,0)R~^(s,0) I = Sup lR(t,0)R*^(s,0)x(s) I 

I x( s) } 

= 1 x( t) I < L 

where L s= C. 

Similarly for s + r < t < s+2r 

lR(t,0) r“^(s,0)! < 1 ? 
and in general for t > s 

lR(t,0) Rr^(s,0)| < e^ < L e^ 

If we take r = ^ / then for t > s+r^ 

L < r(t-s) 

and therefore iR(t,0) R ^(s/0)i < L e^^^ • 

Conversely, ^en (5.4.6) is true, we can have bounded growth 

yx* 

of the syst«n by taking C = max {1,L e } 
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RCTiark 5.4.1 

Lonma 5.4.2 shows that the condition for bounded 
growth is independent of the choice of r. 

Remark 5.4.2. 

The conditions (5*4.5) for exponential dichotomy can 
be expressed in the following equivalent form. 

|R(t,0) Pgi < M R(s/0)P§1 for t > s, 

(5.4.7) 

lR(t,0)(I-P)|l < N e‘■^^^-®^R(s,0)(I-P)§l for s > t, 

and lR(t,0) PR''^(t,0)i < M for all t > O. 

Lemma 5.4.3. 


If (5.4.1) has exponential dichotomy on a subinterval 
then it also has exponential dichotomy on the half line 
iR"**, with the same projection P and the same exponents a and 0. 


Proof : Choosing W = Sup IRCtC/O)! 

0<T<to 

we have lR(t,0) r“^(s.O) i < W for O < s, t < t^. 


If 0 < s < tQ 5 t, then 


|R(t,0)PR'"^(s,0)l = lR(t, 0 )PR“'^(t^, 0 )R(tQ, 0 )R -"(5^0)1 


.- 1 , 


= WtR(t,0) FR” (tQ/0)l 
< WLe ^ ® 


If O < s < t < t^ then 
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lR(t,0)PR“^(s,0)} = |R(t,0)R“^(tQ,0)R(tQ,0)PR“^(tQ,0) 

RCt^^O) r“^(s,o)i 

< lR(tQ,0) PR-^Ct^^O)! 

< W^M < W^H e ° g--a(t-s) ^ 

Hence ! R(t,0)PR“^(s,0) | < M for O < s < t < « 

where 

fo CL't^ „ 

M = max iVMe , VTM e ) . 

Similarly, lR(t,0)(I-p) R~^(s,0)i < L forO<t<s < « 

Now we shall prove a result which indicates that 
bounded growth leads to exponential dichotomy. 

Theorem 5.4.1. 

Suppose that every solution x(t) of (5.4.1) satisfies 
the following conditions 


I x(t) I < Cl x(s) 1 for 0 5. s 5. t < s+T 

(5.4.8) 

and I x( t) t < e Sup | x(u) | for t > T 
“ I u-tl < T “ 

where T, C,e are constants such that T > O, C > 1 and O < 8 < 1 

4. 

Then ec^ation (5.4,1) has exponential dichotomy on 3R . 

Proof : Supposing that x(t) is a non- trivial bounded solution 
then define for s > O, 


X(s) = Sup 1 x(u) 1 
u>s 
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Then for t > s+T 

ix(t)i < e Sup ix{u)i < e x(s) 
lu-tl<T "" 

and therefore Kis) = Sup i x(u) f , 

s<u<s+T 

Further, it follows that 

I x( t) 1 < C 1 x( s) I for 0 < s < t < oo. 

If s+nT < t < s+(n+l)T, then 

lx(t)l < Sup lx(u)i <e^Cix(s)i 

1 u-ti <nT “ 

t~s 

< e“^ c e ^ ix(s)i. 

Thus lx(t)| < M e“‘^^'^"®^lx(s) I for 0 < s < t < » 

where M = S > i ^nd a = - > o. 

Now suppose that x(t) is an unbounded solution with 
ix(0)i =1, Then there exist a sequence Ct^^}, "tn ^ O such 
that 

1 x( tj^) 1 = ^ / 1 x( t) I < % for O < t < tjj . 

0 0 

Then T<tj^<t 2 < ••« and tj^ -• «> as n -• » and moreover 
since lx(t )t < 0 Sup ix(u)| and 

0<u<tjj+T 

I x(u) 1 < 0 “^ Jx(t^)l for 0 < u < 

Suppose t < s and < t < 1 s 

(1 < m < n), then 
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lx(t)l < C = 6^“"^ 

< c ix(s)i 

(s-t) 

< c e**^ e ^ ix(s)i . 

Thus lx(t)i < M x(s) I for 1 t < s < «>^ 

where a = - ~ ~ ~ and M = ^, 

Now let be the subspace of the vector space IR^ 
consisting of the initial values of all bounded solutions of 
(5.4*1),^ and let V 2 be subspace ccxnplimentary to Vj^ • For any 
vector I e V 2 such that i^i - 1 , let x(t) = x(t,§) be the 
solution v^ich takes the value g at t = O* Then x(t,^) is 
unbounded and therefore there exists a least value t^ = t^Cg) 

c 

such that lx(tj^,5)| = g' • We claim that the values tj^Cg) are 
bounded* Suppose not then there exists a sequence of unit 
vectors e V 2 such that -* «» as v ■♦ <». By the 

compactness of the unit sphere in V 2 we may suppose that 
gy -* g , where || 1 =1* Then 

x(t,§y) - x(t/g) for every t > O. 

Since |x(t/gy)l < ^ for 0 < t < it follows that 

lx(t,§)l < ^ for 0 < t < «> 

which is a contradiction to the fact that § ® ^2* 
there exists > O such that tj^ (g ) < Tj^ for all g #■ and 
therefore every solution xCt) with x(0) e V 2 satisfies 
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lx(t)| < M e |xCs)i for % 5, t < s < «>, 

Conse(5u.ently the first two conditions of (5»4-«7) are satisfied 
on the interval [q^,oo) and third condition is obvious* 
Therefore (5.4.1) has an exponential dichotomy on [ 3 ^, 00 ), 

Thus in view of lemma 5.4*3 the ej^onential dichotcxny of (5 •4*’.) 
is guaranteed onlR**, 

Remark 5*4*3* 

For ordinary differential equations a result similar 
to theorem 5*4*1. has been proved by Massera and 
Schaffer [32]* ' 
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